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Introduction and Motivation 
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Motivation 


Path integral with complex weights appear in many important 
physics: 


e Finite-density lattice QCD, 
spin-imbalanced nonrelativistic fermions 


ə Gauge theories with topological 6 terms 


ə Real-time quantum mechanics 


Oscillatory nature hides many important properties of partition 
functions. 
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Example: Airy integral 
Let's consider a one-dimensional oscillatory integration: 


a= | i T as 
p 3 - 


RHS is well defined only if Ima = 0, though Ai(z) is holomorphic. 
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Figure: Real parts of integrands for a = 1 (x10) & a = 0.5i 
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Contents 


ə How can we circumvent such oscillatory integrations? 
=> Lefschetz-thimble integrations 
[Witten, arXiv:1001.2933, 1009.6032] 


ə Applications of this new technique for path integrals 


> General theorem ensuring that Z € R. 
> Application of the theorem to the SU(3) matrix model 


[YT, Nishimura, Kashiwa, in preparation] 


Yuya Tanizaki (University of Tokyo, RIKEN) Lefschetz-thimble path integral Mar. 18, 2015 @ Stony Brook 


5/31 


Lefschetz-thimble integrations 


Introduction to Lefschetz-thimble integrations 
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Again Airy integral 


Airy integral: 


r3 
@)= [ Fei (F tar), 


The integrand is holomorphic w.r.t x 
=> The contour can be deformed continuously without changing the 
value of the integration! 
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Lefschetz-thimble integrations 


Steepest descent contours 
What is the most appropriate contour for our purpose? 


Re[S(x,a)] should be made as small as possible. 
<> The contour should be perpendicular to Re[S(, a)| = cosnt. 


Steepest descent ones J must satisfy Im[S(x, a)] = const. 
because of the holomorphy. 


m 


if 


with a = exp +i0.1. 


Figure: Contour plots for Re[S(x, a)| 
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Rewrite the Airy integral 
There exists two Lefschetz thimbles J, (o = 1,2) for the Airy 


integral: 
d 3 
Ai(a) = > No [ Pi (3 taz) 


Ng: intersection number of the steepest ascent contour Ke and R. 
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Figure: Lefschetz thimbles 7 and duals K (a = exp(0.1i), exp(mi)) 
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Airy integral & Airy equation 


Let us consider the “equation of motion”: 


dz d ei(2?/3+a2) 


= 0. 
Qn a 


This is nothing but the Airy equation: 


Two possible integration contours 
= Two linearly independent solutions of eom. 
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Lefschetz-thimble integrations 


Generalization to multiple integrals 


Model integral: 


z= | dzı -- - d£n exp S(x). 


What properties are required for Lefschetz thimbles J? 
@ J should be a n-dimensional object in C”. 
© Im[S] should be constant on J. 


Yuya Tanizaki (University of Tokyo, RIKEN) Lefschetz-thimble path integral Mar. 18, 2015 @ Stony Brook 


11/31 


Lefschetz-thimble integrations 


Short note on technical aspects 


Complexified variables (a = 1,..., n): Za = Za + ipa. 


Regard x, as coordinates and p, as momenta, 
so that Poisson bracket is given by 


N a 


a=1,2 


O1 Opa Oxa Opal 
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Short note on technical aspects 


Hamilton equation with the Hamiltonian H = Im[S(z,)]: 


df(z, p) 


a ESI (= -()) 


This is Morse’s flow equation: 


d 


= We can find n good directions for J around saddle points! 
[Witten, 2010] 
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Multiple integral: Lefschetz-thimble method 


Oscillatory integrals with many variables can be evaluated using the 
“steepest descent” cycles Js: 


d"r & = ne | d”z eS), 
h. ee), 


Js are called Lefschetz thimbles, and Im[S] is constant on it. 


Ne. intersection numbers of duals X, and R”. 
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Sign problem in MFA 


Sign problem in MFA and its solution: Formal discussion 
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Sign problem in MFA Formal discussion 


Airy integral with real a 


If a € R, the following partition function is real: 


1 r3 
Zairyla) = z [ ropi (= + ar) : 


However, the integrand is complex-valued. 


Why does this happen? 
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Sign problem in MFA Formal discussion 


Airy integral: Antilinear symmetry 


The system have an antilinear symmetry: 


By summing up in pairs +2, the result becomes manifestly real! 
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Airy integral: Lefschetz thimble 


Extend the linear map x ++ —x to the antilinear map z +> —Z. 


The Lefschetz thimbles form invariants under the antiliner map: 


N 
Figure: efschetz thimble I and duals K (a = exp(0.1i), exp(ri)) 


1) 


3 


| 
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Sign problem in MFA Formal discussion 


Generalization: Set up 


Consider the oscillatory multiple integration, 


z= | d”zx exp — S(x). 


To ensure Z € R, suppose the existence of a linear map L, satisfying 


@ S(x) = S(L- x). 
ə LP =1, 


Let’s construct a systematic computational scheme of Z with 
ZER. 
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Sign problem in MFA Formal discussion 


Generalization: Flow eq. and Complex conjugation 


Morse’s downward flow: 


Complex conjugation of the flow: 


dt Oz; 


and thus z’ := L -Z satisfies Morse’s flow equation . 
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Sign problem in MFA Formal discussion 


Generalization: Immediate conclusion 


Denote the antilinear map as K(z) := L-Z. 
@ When z° is a saddle point, so is K(z7). 
@ Let J, is the Lefschetz thimble around 27. 


=> JE :={K(z) | z € Je} coincides with the Lefschetz thimble 
around K(z”) (up to orientation). 


© Intersection numbers are the same: (Ko, R”} = (KE, R”) under 
appropriate choice of the orientation. 


@ The contributing thimbles always form an invariant pair 


ge. 
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Sign problem in MFA Formal discussion 


General Theorem 


Let us decompose the set of saddle points into three parts, 
X= Yo UN} UX, where 


Xo = {o| = L-27}, 
{o | ImS(27) > 0}, 
x. = {o |ImS(2z7)< 0}. 


M 
+ 
| 


The antilinear map gives X, ~ %_. 
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General Theorem 


The partition function: 


Z = ne | d”z exp —S(z) 
oEX0 To 
+ n f d”z exp —S(z) 
fen, TIA 


Each term on the r.h.s. is real. 


(YT, Nishimura, Kashiwa, in preparation) 
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Sign problem in MFA 


Sign problem in MFA and its solution: 
Application to QCD-like theories 
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penises 2 
The QCD partition function at finite density 


The QCD partition function: 
Zacp = [ra detM (pak, A) exp —Sym [A], 


w./ the Yang-Mills action Sym = $tr m dzt f dæ|F l? (> 0), and 


detM(u, A) = det |y” (8, + ig AL) + Vlak + Max] - 


is the quark determinant. 
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SW HELEN TE Application to QCD 


Charge conjugation 


If Hax # 0, the quark determinant takes complex values 
=> Sign problem of QCD. 


However, Zoco € R is ensured thanks to the charge conjugation 
Ar —At: 


det M(jigk, A) = det M(—jigx, A") 
det M (fax, — A). 


(First equality: y5-transformation, 
Second equality: charge conjugation) 
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SmceumueVigame Application to QCD 


Polyakov-loop effective model 


The Polyakov line L: 


1 re . 
L= ;diag fe 91402), ei -9+92), 02102], 
Let us consider the SU (3) matrix model: 
Zacp = [nano 05) exp [—Ver(41, 65)| ; 


where H = sin? 0, sin? (04 + 305) /2sin? (01 — 305) /2. 
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Smee Vigame Application to QCD 


Charge conjugation in the Polyakov loop model 


Charge conjugation acts as L + Li: 
Vest (2 22) = Verla, 23). 
Simple model for dense quarks (£ := trL): 


(=) 


Va = =f (cl'ls(01,02) + ©*13(1,0,)) 


SY cosh ju (2 cos 6; cos da + cos 202) 


+isinh u (2 cos 6; sin da — sin 20,) k 
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SWEET Application to QCD 


Behaviors of the flow 


Olt pe eee eee ee He ee 


[6] 2 


J 
Re(z1) 
Figure: Flow at h = 0.1 and u = 2 in the Re(z1)-Im(z2) plane. 
The black blob: a saddle point. 
The red solid line: Lefschetz thimble 7. 
The green dashed line: its dual K. (YT, Nishimura, Kashiwa, in preparation) 
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SruwelelJlumlm Wi Application to QCD 


Saddle point approximation at finite density 


The saddle point approximation can now be performed. 


Polyakov-loop phases (z1, z2) takes complex values, so that 


(0), ® ER. 


Since Im(z2) < 0 and 
£ > —(2e'” cos 6, + e 7”), 


we can confirm that 


(YT, Nishimura, Kashiwa, in preparation) 
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Summary 


ə Lefschetz-thimble integral is a useful tool to treat multiple 
integrals. 


ə Saddle point approximation can be applied without violating 
ZER. 


ə Sign problem of effective models of QCD is (partly) explored. 


Yuya Tanizaki (University of Tokyo, RIKEN) Lefschetz-thimble path integral Mar. 18, 2015 @ Stony Brook Bigs 


